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SUMMARY 


The  nns  error  in  the  source  bearing  estimates  was  observed  to  in- 
crease with  decreasing  signal-to-noise  ratio  (SNR),  when  the  model 
order  (p)  approached  the  number  of  array  elements  (N  = 31) , and  when 
the  source  of  interest  was  not  well  separated  from  adjacent  sources. 
Typical  rms  errors  observed  for  isolated  sources  (for  p - N/3)  ranged 
from  0.02  A/D  at  a signal-to-noise  ratio  (SNR)  of  17  db  to  0.5  A/D  at 
-9  dB.  The  statistical  variability  and  number  of  spurious  peaks  in  the 
array  response  were  observed  to  increase  rapidly  as  p approached  N.  At 
an  SNR  of  less  than  0 dB  the  array  response,  even  when  averaged  over  50 
independent  data  sets,  indicated  numerous  spurious  sources  when  the 
model  order  was  large.  The  equivalent  number  of  degrees  of  freedom  (df) 
of  the  array  response,  estimated  as  twice  the  square  of  the  average 
divided  by  the  variance  of  the  response,  were  calculated  and  compared 
with  df  ■ N/p  as  conjectured  by  Parzen  (1970)  and  demonstrated  (asymp- 
totically) by  Kromer  (1969).  Parzen' s conjecture  was  consistent  with 
the  observed  results  for  p <<  N,  but  was  found  to  be  very  optimistic 
at  high  model  orders.  It  was  also  found  that  the  spatial  response  to  an 
individual  point  source  increases  in  width  in  the  presence  of  nearby  tar- 
gets as  well  as  with  decreasing  signal-to-noise  ratio;  and  while  the 
spatial  response  to  a single  point  source  is  typically  much  narrower  for 
maximum  entropy  processing,  the  minimum  angular  separation  required  for 
the  resolution  of  adjacent  sources  is  only  slightly  less  (roughly  a fac- 
tor of  one-half)  than  that  of  a conventional  beamformer  even  at  high 
signal-to-noise  ratios. 
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MOTIVATION 

B<„.iC°n8Id!r  thC  problem  of  estimating  the  bearing  and  strength  of  a 
single  point  source  using  a linear  array  where  the  source  may  either 
be  actively  generating  Its  own  signal  or  reflecting  a transmitted  wave- 
Assuming  the  source  and  receiver  are  separated  sufficiently  to 
approximating  the  received  signal  as  a plane  wave,  the  fields 
in  the  vicinity  of  the  array  may  be  written  as 

s(x,t)  - p (t-a*x)  (D 

where  a - v/|v|2,  v is  the  velocity  vector  of  the  plane  wave,  and  x 
specifies  the  field  point.  Further,  if  the  signal  is  assumed  to  be 
narrowband,  then 


_ -i(u>  a*x-u>  t) 

8(x,t)  - w(t-a*x)e  0 (2) 

_ _ -io)  a»x 

***!;*  °-  18  the  comPlex  envelope  of  the  waveform  observed 

at  the  field  point  x.  The  array  geometry  is  shown  in  Figure  1.  Assuming 
the  linear  array  consists  of  N elements  spaced  Ax  apart  and  lies  along 
the  x axis  of  the  coordinate  system  with  the  origin  at  one  end,  then  the 
output  of  the  nth  element  is 


Zn(t)  " s(nAx  x,t) 


w(t  - nAx  a )e 
x 


-12  irk  nAx 


e'V  (3) 


where  x is  a unit  vector  along  the  x-axls. 


kx  " Vx  ‘ I 8in  (0> 

Is  the  X component  of  the  vector  wave  number,  6 is  the  angle  between  the 
direction  of  propagation,  and  the  x axis,  X Is  the  wavelength,  and  a - 
a*x  " a *8  the  x component  of  the  Inverse  velocity  vector.  x 

Equation  (3)  may  be  approximated  as 


*n(t)  91  w(t)  e 


-i2irk  nAx  iu  t 
X e ° 


(Text  Continued  on  Page  3) 
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It  the  transit  time  .uross  the  arrav  is  small  compared  to  the  time  of 
variation  of  the  signal  envelope;  i.e.,  tin*  inverse  of  the  signal  band- 
width. Finally,  the  carrier  term  in  Equation  (•'*)  nuv  be  removed  bv 


coherently  detecting  z (t)  to  obtain 
‘ n 


z’  (t) 
n 


W(  t h 


• id  k nAx 
x 


(5) 


All  of  the  information  relevant  to  the  source  strength  and  bearing 
is  contained  in  the  z’(t)  and  their  relative  delay  times*'1'.  The  above 
equation  Indicates  that,  under  the  assumptions  made,  the  fields  measured 
along  a linear  array  at  time  t vary  sinusoidally  as  a function  of  x with 
frequency  kx  cycles/unit  distance.  Knowledge  of  k^  and  w , are  sufficient 
to  determine  the  direction  of  propagation  in  a noiul isperssi ve  medium  when 
the  velocity  of  all  possible  plane  waves  are  required  to  lie  in  the  xy 
plane  since  may  In'  obtained  from  the  dispersion  relation 


(2n)-  ,ko 


k + k 

X V 


(tO 


where  c is  the  velocity  of  propagation.  Thus,  estimating  the  bearing  to 
a single  point  target  is  nothing  more  than  estimation  of  the  frequency 
of  the  variation  along  the  array.  It  is  evident  from  the  discussion 
above  that  the  problem  is  simply  one  of  spectral  estimation.  The  conven- 
tional approach  to  this  problem  is  beamforming  which  for  the  narrowband 
case  reduces  to  calculation  of  the  Discrete  Fourier  Transform  of  the 
array  element  outputs  at  a particular  time.  The  search  for  high  spatial 
resolution  quite  naturally  leads  to  algorithms  developed  for  high  reso- 
lution spectrum  analysis. 

In  an  active  sonar  system  the  coherently  detected  array  outputs  are 
a superposition  of  terms  of  the  form  given  in  liquation  (5),  whore  w(t) 
is  an  appropriately  delayed  version  of  tin*  transmitted  signal  envelope. 
The  principal  problem  is  to  resolve  scattering  centers  lying  at  roughly 
the  same  range  (more  precisely  in  the  same  range  cell  where  a range  cell 
is  one-half  the  velocity  of  propagation  times  the  Inverse  of  the  trans- 
mitted signal  bandwidth)  since  scattering  centers  lying  at  differing 
ranges  mav  be  distinguished  in  the  time  domain. 

In  a passive  system  the  problem  is  to  resolve  sources  emitting 
signals  which  overlap  in  time.  It  is  important  to  improve  the 


^*'1.  N.  El-Beherv  and  R,  H.  MacFhie,  .'fax  tmum- 1 i kol  ihood  Estimation  of 
Source  I'a ramct  ers  from  Time  Sampled  Outputs  of  a linear  Arrau,  d. 
Acoust.  Soc . Am.  b2,  p.  125-114  (1977). 
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signal-to-noise  (SNR)  as  much  as  possible  by  coherent  integration  prior 
to  spectral  (bearing)  estimation  as  t tie  maximum  entropy  method  is  quite 
sensitive  to  noise.  In  either  ease  the  appropriate  coherent  integration 
times  are  limited  in  t lie  act ive  ease  by  transmitted  signal  bandwidth,  and 
in  the  passive  ease  by  the  frequency  stability  of  the  Source.  In  many 
cases  it  is  possible  to  improve  the  estimates  still  further  bv  averaging 
the  results  ot  multiple  independent  observations  obtained  bv  appropriately 
sampling  the  coherently  detected  array  outputs.  This  should  be  used 
whenever  possible. 

Ihe  principal  interest  was  in  active  high  resolution  sonar  systems 
for  the  purpose  of  localizing  and  resolving  discrete  targets  where  only 
a single  observation  of  the  coherently  detected  array  outputs  was  avail- 
able. In  this  report  maximum  entropy  processing  as  a means  of  locating 
and  resolving  point  sources  based  on  a single  time  sample  (observation) 
of  the  array  outputs  will  be  investigated.  Derails  of  the  maximum 
entropy  method  will  be  discussed  only  briefly  in  the  following  section 
since  these  are  readily  found  in  the  literature.  the  principal  results 
of  this  study  are  presented  in  the  SIMULATION  AND  RKSULTS  section. 

This  section  begins  with  a discussion  of  the  choice  of  model  order,  and 
its  relation  to  the  statistical  stability  of  the  array  response.  This 
is  followed  by  a discussion  of  the  effects  of  SNR  and  model  order  on  (1) 
errors  in  bearing  estimation,  (2)  the  width  of  the  array  response  to 
individual  sources,  and  (1)  the  minimum  separation  required  for  reso- 
lution of  equal  amplitude  point  sources.  When  possible  our  results  are 
compared  with  those  of  other  researchers . The  report  concludes  with  a 
summarv  of  the  principal  findings. 


MAX  I MIDI  ENTROPY  SPECTRAL  ANALYSIS 


The  maximum  entropy  spectral  analysis  technique  was  first  proposed 
by  J.  P.  Burg  in  1967  as  a means  of  achieving  high  resolution  spectral 
estimates^ . This  approach  was  later  recognized  to  be  identical  to 
auto-regressive  spectral  analysis,  linear  predictive  deconvolution 
filtering,  and  all  pole  modeling  of  the  spectrunC'.  The  entropy  of  a 
Gaussian  band  limited  time  series  is  proportional  to 


Burg , J.  P.,  ’fa i mum  Entropy  Spt'ctral  Analysis,  presented  at  the  37th 
Annual  International  Mooting,  Soc.  of  Explor.  Geophvs.,  Oklahoma  Citv, 
Oklahoma  (1967). 

N.  McDonough,  Maximum  Entropy  Spatial  Processing  of  Arrau  Data, 
Geophysics  39,  p.  843-851  (1974). 
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where  Sit)  is  t (to  power  speetrutn,  uni  f is  the  Nvquist  frequency^  . 

The  power  spectrum  is  the  Fourier  transform  of  t Uv-  lutocorre  1 at  ion  func- 
tion ot  tin'  time  series  which  is  known  onlv  over  a finite  interval. 

Burs  suggested  that  the  most  reasonable  choice  for  the  unknown  auto- 
correlations is  the  one  which  maximizes  the  entropy.  Maximization  of  H 
subject  to  the  constraints  imposed  bv  the  known  autoeorrelat ions  im- 
plies .1  power  spectrum  of  the  form 


■i-’-  f i ' 1 1 


where  1,  a*[ , aC.-a,  are  the  coefficients  of  t lie  p+1  length  prediction 
error  filter,  and  p'  is  the  prediction  error'  . Burg  later  suggested 
an  approach  by  which  the  filter  coefficients  could  be  estimated  directly 
from  the  data'  . This  approach  is  based  on  choosing  the  reflection  co- 
etticients  (a!-,)  to  minimi/.c  the  prediction  error  power  obtained  from 
running  the  filter  both  forward  and  backward  across  the  data  and  on  the 
Levinson  recursion  algorithm,  and  is  outlined  in  Reference  b.  In  this 
study  the  maximum  entropv  array  response  (spatial  spectrum)  was  calcu- 
lated using  Equation  (8)  with  the  replacements  f *■  Ax,  and  At  *•  Ak  and 
using  the  coherently  demodulated  element  outputs  (at  a specified  t^me) 
as  the  sampled  complex  data  series.  Burg's  method  (modified  for  use 
with  complex  data)  of  coefficient  calculation  was  used  to  obtain  the 
results  presented  in  this  report  because  it  is  widelv  used,  and  has 
been  reported  to  produce  resolution  superior  to  the  Yule  Walker 


Smvlie,  D.  F.  , Clarke,  0.  K.  C.  and  I’lrvch,  I'.  J.,  .Ana  it/s  is  of 
Irreaular : ties  in  the  earth’s  R otat ion,  in  B.  Alder,  S.  Fernback, 
and  M.  Rotenberg  (Eds.)  "Methods  in  Computational  1'hysies,  Yol.  Id," 
New  York:  Academic  Press,  pp.  191-410  (1971). 

' 1 ' Burg , d.  P.  , A .Vow  Analysis  Technique  for  riru'  Series  Data,  presented 
at  NATO  Advanced  Study  Institute  on  Signal  Processing,  Enschede, 
Netherlands  (1968). 

*’  • Andersen , N.  , On  the  Calculation  of  Filter  Coefficients  for  Maximum 
Entropy  Spectral  Analysis,  Ooophysics,  Vol.  19,  No.  1,  pp.  69-71, 
February  197-., 
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approach  .It  has,  however,  been  pointed  out  by  Makhoul  that  the 

Burg  method  of  calculating  the  reflection  coefficients  is  only  one  of 
a number  of  "suboptimal"  solutions  to  the  lattice  minimization  prob- 
lem^1*.  The  complex  formulation  of  the  maximum  entropy  algorithm  is 
discussed  in  References  3 and  10. 


SIMULATION  AND  RESULTS 


In  the  following  simulations  a linear  array  has  been  assumed  having 
31  elements  spaced  one-half  wavelength  apart  responding  to  inphase  mono- 
chromatic signals  generated  by  two  point  sources.  The  inphase  condition 
was  chosen  because  this  is  the  most  difficult  condition  to  resolve.  The 
signals  observed  by  each  element  are  assumed  to  be  corrupted  bv  narrow- 
band  Gaussian  noise  which  is  independent  between  elements.  After  co- 
herent demodulation  the  output  of  the  nth  element  is 


z’  (t) 
n 


2 -i2TtkCnAx 

Z w"(t)  e x + n (t) 
i n 

«.« I 


(9) 


where  n (t)  is  the  envelope  of  the  complex  noise  added  to  the  nth  ele- 
ment. ¥he  SNR’s  quoted  in  the  remainder  of  this  study  are  defined  as 
the  ratio  of  the  signal  power  (generated  by  one  of  the  sources)  to  noise 
power  at  the  output  of  each  element  after  demodulation. 

At  each  combination  of  SNR  and  target  separation  50  data  sets  (with 
independent  noise  realizations)  corresponding  to  a single  time  sample  of 


<3> ibid. 

<7>T.  J.  Ulrych  and  T.  N.  Bishop,  Maximum  Entropy  Spectral  Analysis  and 
Autoregressive  Decomposition,  Reviews  of  Geophvsics  and  Space  Phvsics 
13,  p.  183-200  (1975). 

R.  Radoski,  P.  F.  Fougere,  and  E.  J.  Zawalick,  A Comparison  of 
Power  Spectral  Estimates  and  Applications  of  the  Maximum  Entropy 
Method,  J.  Geophysical  Research  80,  p.  619-625  (1975). 

J.  Makhoul,  Lattice  Methods  in  Spectral  Estimation,  in  the  "Proceed- 
ings of  the  RADC  Spectrum  Estimation  Workshop,"  p.  159-173  (1978). 

Texas  Instruments,  Inc.,  Technical  Report  No.  1,  Advanced  Signal 
Processing,  ALEX(03)-TR-75-01 , The  Maximum  Entropy  Spectrum  and 
the  Burg  Technique,  by  T.  E.  Barnard,  1975. 
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each  of  the  N array  elements  were  examined.  Coefficients  corresponding 
to  model  orders  1 through  30  were  calculated  for  each  data  set.  Spectra 
were  calculated  corresponding  to  model  orders  of  5,  10,  15,  20,  25,  and 
30.  The  reported  angular  separations  were  achieved  by  holding  the  bear- 
ing of  one  target  fixed  and  moving  the  second  target  relative  to  the 
first.  This  was  done  so  that  changes  in  the  response  to  one  of  the 
targets  would  be  related  to  changes  in  target  separation  and  not 
changes  in  target  bearing.  Only  parameters  of  the  array  response  rela- 
ted to  the  fixed  target  are  reported.  Some  definite  trends  were 
observable  in  the  results. 


CHOICE  OF  MODEL  ORDER 


The  estimation  of  the  order  of  the  appropriate  maximum  entropy 
model  has  been  treated  by  a number  of  authors.  Based  on  the  results 
reported^  ^ by  Ulrych  and  Bishop,  Akaike's  Final  Prediction  Error  (FPE) 
Scheme  was  selected  in  which  the  optimum  value  of  p 
value  which  minimizes  the  FPE  given  by 


is  chosen  as  that 


(FPE)  = ^-2.  p 
P N - p Kp 


Unfortunately,  in  this  application  the  minimum  of  the  FPE  does  not  ap- 
pear to  indicate  acceptable  model  orders.  Typical  results  for  the 
average  and  standard  deviation  of  the  FPE  over  50  data  sets  for  model 
orders  1 through  30  are  shown  in  Figure  2.  Two  characteristics  of  the 
FPE  are  quite  noticeable.  First,  the  variance  of  the  FPE  increases  with 
model  order.  This  would  appear  to  be  related  to  the  fact  that  the  sta- 
tistical stabilitv  of  the  estimated  spectra  decreases  with  increasing 
model  order.  And  second,  the  FPE  has  a local  minima  at  below  10,  is 
relatively  flat  up  to  p = 25,  and  decreases  rapidly  for  p > 25.  Based 
on  the  average  FPE,  p » 30  would  be  chosen,  however  a comparison  of  the 
spectra  obtained  for  differing  model  orders  indicates  that  this  is  not 
acceptable  since  both  the  number  of  spurious  peaks  and  statistical  vari- 
ability of  the  resulting  spectra  increase  rapidly  as  p approaches  N. 
Clearly  if  the  calculated  responses  are  to  be  meaningful  both  the  num- 
ber of  spurious  peaks  and  the  statistical  variability  of  the  response 
must  be  controlled.  This  is  particularly  important  when  the  responses 
are  based  on  single  observations  of  the  array  outputs  and  no  averaging 
is  performed.  These  problems  are  more  acute  at  low  SNR  because  the 
signal  peaks  are  relatively  smaller,  but  they  can  be  significant  even  at 
high  SNR  when  p is  large.  Even  the  average  spectrum  may  be  of  marginal 
value  when  the  model  order  approaches  N.  Figures  3a  and  3b  compare  the 
average  (over  50  data  sets)  spectrum  at  an  SNR  of  7 dB  for  p = 30,  and 
p » 15;  clearly  the  spectrum  is  useless  for  p = 30  due  to  the  number  of 
spurious  peaks.  At  an  SNR  of  -3  dB  the  largest  model  order  producing  a 


(7)ibid. 
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Model  Order  15 


Bearing  (sin  0) 


A 

o 


FIGURE  3a.  AVERAGE  ARRAY  RESPONSE 


useful  average  wt.a  p - JO.  In  «nv  event,  i(  appears  that  I It  tin  Is 
gained  In  terms  at'  resolution  or  hearing  estimation  aeeuraev  hv  choosing 

p ' 10. 


As  In  or Jinary  sport ra!  estimations,  t In*  equivalent  number  of 
degrees  oi  freedom  of  the  arrav  response  may  be  estimated  as 


v 


\ Jir(s(kx))| 

j V-'7rTs(k*T)  j 


over  all  k 

x 


(10) 


where  K(8(kx))  and  Var(S(kx))  are  the  mean  and  varlanee  of  the  maxi- 
mum entropy  array  response.  This  quantity  provides  a measure  of  the 
statistical  stability  of  the  response.  Figure  s Illustrates  the  depen- 
dence of  this  quantity  on  the  model  order  where  the  degrees  of  freedom 
were  estimated  from  the  mean  and  varlanee  ot  the  response  over  the  r»0 
data  seta  for  the  SNK  and  target  separation  Indicated,  This  behavior 
is  typical  of  the  results  observed.  It  Is  apparent  that  the  statis- 
tical variability  of  the  array  response  Increases  rapidly  as  p ap- 
proaches N.  There  are  few  theoretical  results  on  the  statistical 
properties  of  autoregressive  spectral  estimates  for  comparison  with 
these  results.  it  has,  however,  been  conjectured  bv  Parzen  (1970)^11* 
and  demonstrated  bv  Kromer  (1969)^1“^  that  aut  ogregresslve  modeled 
spectral  estimates  tire  asvmpf ot 1 ca I I y (for  large  p and  N with  p N) 
distributed  with  v ■*  N/p  degrees  of  freedom^13'.  The  statistical  sta- 
bility of  the  observed  results  Is  consistent  with  Parzen's  conjecture 
when  p *•<  N,  but  Is  increasingly  less  stable  as  p approaches  N as  was 
expected  from  the  asymptotic  nature  of  the  theoretical  result. 


It  was  concluded  from  observations  that  the  best  combination  of 
statistical  stability  and  source  resolution  was  obtained  for  p 1 10 
although  p - 5 inav  be  acceptable  It"  the  resolution  requirements  are  not 
stringent.  Choosing  p - 1/3  the  number  of  array  elements  appears  to  be 
a good  starting  point  although  the  "optimum"  choice  depends  on  the  SNR, 


^ll*Parzon,  K. , Multivariate  Time  Series  Modelin>i,  in  "Multivariate 

Time  Series  Analysis  II,"  K.  Kirshnlah.  lid.,  Academic  I’ress,  New 
York,  1970. 

^ 1 Kromer , R.  li. , Asymptotic  I'roperties  of  the  Autoregressive  Spectral 
Estimator,  Olssertat ton.  Dept,  of  Statistics,  Stanford  University 
1969. 

^ 1 •**  llersch,  W.  and  l.lu,  R.  11.  , Autoregressive  Model  Spectral  Estima- 
tion, Some  Simulation  Studu  Statistical  I'erformance  Results,  in 
"Proceedings  of  the  RADC  Spectrum  list l mat  l on  Workshop,"  May  1978. 
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number  of  data  records  available,  and  the  desired  resolution,  as  well 
as  the  nature  of  the  data.  It  is  interesting  to  note  that  the  choice 
of  p - 10  gives  roughly  the  same  stability  as  that  obtained  with  con- 
ventional beamforming  with  an  unweighted  aperture  (conventional  spec- 
tral estimation  with  rectangular  window). 


BEARING  ERROR 

Source  localization  in  azimuth  is  particularly  important  in  many 
applications.  It  appears  that  the  maximum  entropy  procedure  provides 
excellent  bearing  estimates  when  the  source  of  interest  is  well  sepa- 
rated from  other  sources. 

Source  bearing  was  estimated  by  locating  the  largest  peak  in  the 
vicinity  of  the  known  source  location.  Table  1 indicates  the  total  rms 


TABLE  1 


RMS  DEVIATION  OF  PEAK  LOCATION  FROM  TRUE  BEARING 


RHB 

Separation* 

RMS  Deviation* 

Model  Order 

1 - 5 1 

1 IQ 1 

1 is  1 

-9 

CO*  * 

mm 

mm 

m 

1.5 

BA 

ESI 

2.0 

■sHU 

0.9 

0.6 

-3 

CO 

0.2 

0.2 

BB 

1.5 

0.8 

0.6 

i EK# 

2.0 

1.2 

1 

0.2 

i 

0.2 

7 

CD 

mam 

ESI 

1.5 

111 

ESI 

2.0 

la 

■i 

■a 

17 

CO 

0.02 

ESI 

1.5 

0.2 

mxm 

2.0 

0.05 

Kb 

KB 

*in  multiples  of  A/D 
**single  source 
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deviation  of  the  estimated  bearings  from  the  true  bearing.  The  total 
deviation  error  is  a function  of  the  model  order,  SNR,  and  target  sepa- 
ration. It  tends  to  increase  with  decreasing  SNR  when  the  model  order 
is  large  and  when  the  targets  are  not  well  resolved. 

The  total  deviation  error  may  be  interpreted  as  consisting  of  two 
parts:  the  bias  or  average  error  (over  the  50  data  sets)  in  source 

location  and  deviations  around  this  average.  The  bias  is  not  signifi- 
cant when:  only  a single  source  is  present,  appears  to  be  relatively 
unimportant  for  large  separations,  and  tends  to  increase  as  the  separa- 
tion decreases.  At  high  SNR  the  bias  is  the  major  source  of  error  when 
multiple  sources  are  present.  The  deviations  around  the  average  increase 
as  the  SNR  decreases  and  is  the  major  source  of  error  for  SNR  < 0 dB 
(provided  the  sources  are  resolved).  It  is  interesting  to  note  from 
Table  1 that  for  isolated  sources  the  errors  for  p = 5 are  not  signifi- 
cantly greater  than  for  p = 10  or  15.  If  resolution  of  multiple  sources 
is  not  required,  p = 5 may  be  acceptable  for  source  localization,  and 
as  noted  previously  gives  a considerably  more  stable  response. 


It  is  instructive  to  compare  the  isolated  target  results  with  the 
Cramer-Rao  lower  bound  on  the  rms  error  given  (in  multiples  of  X/D) 
by 


ii) 


1 /1+N('SNR)  /6 

O » — * 

vQ  SNR  2tt,^ 


where  Q is  the  number  of  independent  time  samples  used  in  the  estimate 
(in  this  case  Q = 1).  The  C-R  lower  bound  is  the  limit  for  the  per- 
formance of  any  estimator.  For  SNR's  of  17,  7,  -3,  and  -9  dB  the  result- 
ing C-R  bounds  are  1x10“%  3x10"%  1x10"^,  and  2xl0“l.  The  results  in 
Table  1 for  p « 10  are  roughly  twice  the  C-R  bounds.  This  appears  to 
compare  reasonably  well  with  the  results  on  maximum  liklihood  estimation 
reported  in  Reference  1. 


PEAK  WIDTHS 

Provided  the  targets  were  resolved,  the  3 dB  width  of  the  response 
to  the  fixed  source  was  always  much  less  than  the  response  width  for  con- 
ventional beamforming  of  0.9  X/D.  Even  at  SNR  = 0 dB  and  p = 5 the 
average  response  width  to  a single  source  was  only  0.5  X/D.  The  width 
increased  as  the  SNR  or  model  order  decreased,  and  as  target  separation 
decreased. 


(1)  ibid. 
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It  is  interesting  to  compare  the  observed  single  target  response 
widths  with  those  which  can  be  calculated  theoretically.  The  expres- 
sion for  the  3 dB  response  width  for  maximum  entropy  array  processing 
may  be  obtained  directly  from  the  results  obtained  by  Lacoss^14^  in  the 
time  domain  if  the  analogy  between  k and  f is  noted,  and  between  Ax 
and  At  x 

Akx  “ irAxp2  (SNR)  (11) 

where  kx  is  the  3-dB  (wave  number)  width  of  the  spatial  response.  This 
is  easily  related  to  the  angular  width  A6  since 

k * 7-  sin  6 
x A 

Ak  » y cos  A0  (12) 

X A 

2N  A 

irp^(SNR)  cos©  V 

Table  2 compares  the  predicted  and  observed  angular  widths  for  the 
three  SNR' s when  0=0  (source  located  on  the  perpendicular  bisector  of 
the  array  axis).  Substitution  of  the  appropriate  model  order  predicts 
widths  somewhat  less  than  those  actually  observed.  This  is  not  surpris- 
ing because  Equation  (11)  was  derived  assuming  exact  knowledge  of  the 
signal  autocorrelations  rather  than  estimates. 


TABLE  2 

BEAMWIDTH  COMPARISON 
Observed* 


SNR  (dB) 

p - 5 

p ■ 10 

p * 15 

17 

8 x 10~3 

4 x 10-3 

3 x 10'3 

7 

CM 

1 

O 

H 

X 

3 x 10-2 

2 x 10~2 

-3 

9 x 10_1 

3 x 10"1 

5 x 10'2 

-9 

o 

o 

H 

X 

CM 

1 x 10° 

5 x 10"1 

Calculated 

3.9  x 10_1/p2 

0 2 

3.9  x 10  /p 

3.9  x lO1/?2 

1.6  x 102/p2 


*in  multiples  of  A/D 


(14)Lacos8,  R.  T. , Data  Adaptive  Spectral  Analysis  Methods,  Geophysics, 
Vol.  36,  No.  4,  pp.  661-675  (August  1974). 
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It  should  be  noted  that  while  the  resolution  of  a conventional 
beamformer  is  directly  related  to  its  response  to  a single  point  source, 
this  is  not  the  case  for  maximum  entropy  processing  since  the  algorithm 
is  nonlinear  and  superposition  does  not  hold.  It  is  sometimes  stated 
that  the  maximum  entropy  spectrum  has  no  sidelobes.  This  can  be  mis- 
leading since  it  applies  only  to  an  isolated  source.  There  is,  in  fact, 
a nonlinear  interaction  between  components  as  evidenced  by  the  changes 
in  response  width  as  a function  of  target  separation  which  leads  to 
problems  similar  to  those  caused  by  sidelobes  in  conventional  beamform- 
ing. Based  on  the  observed  single  target  response  it  might  be  concluded 
that  the  "resolution"  was  many  orders  of  magnitude  better  than  is  actu- 
ally the  case. 

Resolution  is  a difficult  parameter  to  define  for  maximum  entropy 
spectra.  When  the  target  separation  is  small,  it  is  often  difficult  to 
determine  if  the  observed  peaks  correspond  to  sources  or  are  due  to 
noise.  For  the  purposes  of  this  study  equal  amplitude  targets  were  con- 
sidered to  be  resolved  when  there  was  a dip  in  average  (over  50  samples) 
response  of  3 dB  between  peak  responses  to  the  targets.  In  general,  the 
results  were  as  follows:  the  minimum  resolvable  separation  (MRS)  in- 
creased as  the  SNR  or  model  order  decreased.  The  decrease  in  MRS  with 
SNR  may  be  partially  compensated  for  by  increasing  the  model  order, 
however  as  noted  previously,  the  number  of  spurious  spectral  peaks  in- 
creased rapidly  as  the  model  order  increased.  For  model  order  10  the 
minimum  resolvable  separation  (MRS)  for  SNR's  of  17,  7,  -3,  and  -9  dB 
were  0.7,  1.0,  1.5,  and  2 A/D.  The  MRS  for  model  order  15  was  only 
slightly  better  while  the  MRS  for  a model  order  of  5 was  roughly  twice 
that  obtained  with  a model  order  of  10.  In  cases  where  the  resolution 
requirements  are  not  extreme  a low  model  order  may  be  acceptable. 

For  comparison  the  minimum  angular  separation  required  at  a high 
SNR  for  the  resolution  of  two  equal  amplitude  in  phase  point  target* 
using  conventional  (delay  and  sum)  beamforming  if  1.4  A/D,  An  exten- 
sive study  of  the  resolution  of  sinusoids  is  described  by  Marple 
(1976) ^1B'.  The  average  MRS  reported  is  roughly  one-half  that  reported 
above^16).  While  there  are  differences  in  the  signals  (real  versus  com- 
plex) and  the  methods  employed,  it  appears  that  this  discrepancy  is 
primarily  caused  by  the  differing  definition  of  resolution  used  in  the 
current  study. 


Marple  Conventional  Fourier,  Autoregressive,  and  Special  ARM 
Methods  of  Spectrum  Analysis,  Engineers  Dissertation,  Department 
of  Electrical  Engineering,  Stanford  University  (1976). 
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